arXiv:1508.05379vl [hep-ex] 21 Aug 2015 



ELSEVIER 


Available online at www.sciencedirect.com 


ScienceDirect 


Nuclear Physics B Proceedings Supplement 00 (2015) l-[^ 


Nuclear Physics B 
Proceedings 
Supplement 


Low-energy (anti)neutrino physics with Borexino: 
Neutrinos from the primary proton-proton fusion process in the Sun 


P. Mosteiro'* ', G. Bellini‘S, J. Benziger'', D. Bick", G. Bonfini^, D. Bravo®, B. Caccianiga'^, L. CadonatE’^, 

F. Calaprice^’*, A. Caminata^, P. Cavalcante^, A. Chavarria^, A. Chepurnov^, D. D’Angelo‘S, S. Davini% A. Derbin”, 
A. EmpP, A. Etenko'‘’j, K. Eomenko^’^, D. Eranco*^, E. Gabriele^, C. Galbiati'', S. Gazzana^, C. Ghiano^, 

M. Giammarchi'^, M. G6ger-Neff°, A. Goretti^*, M. Gromov**, C. Hagner", E. Hungerford^ Al. lanni**. An. lanni**, 
V. Kobychev*, D. Korablev'*, G. Korga^ D. Kryn**, M. Laubenstein^, B. Lehnert*, T. Lewke**, E. Litvinovich’^J, 

E. Lombardi**, P. Lombardi**, L. Ludhova**, G. Lukyanchenko’^, I. Machulin’^A S. Manecki®, W. ManeschgP, 

S. Marcocci*’", Q. MeindP, E. Meroni**, M. Meyer", L. Miramonti", M. Misiaszek', M. Montuschi**, V. Muratova", 
L. Oberauer**, M. Obolensky**, E. Ortica*<, K. Otis**, M. Pallavicini*, L. Papp®’**, L. Perasso*, A. Pocar**, G. Ranucci", 
A. Razeto**, A. Re**, A. Romani*!, jsj Rossi**, R. Saldanha**, C. Salvo*, S. Schonert**, H. SimgenP, M. Skorokhvatov**’!, 
O. Smirnov**, A. Sotnikov**, S. Sukhotin**, Y. Suvorov*"’****, R. Tartaglia**, G. Testera*, D. Vignaud**, R. B. Vogelaar®, 
E. von Eeilitzsch**, H. Wang"*, J. Winter**, M. Wojcik*, A. Wright**, M. Wurm**, O. Zaimidoroga**, S. Zavatarelli*, 

K. Zuber*, G. Zuzel* 

^Physics Department, Princeton University, Princeton, NJ 08544, USA 
^ Joint Institute for Nuclear Research, Dubna 141980, Russia 
^Dipartimento di Fisica, Universitd degli Studi e INFN, Milano 20133, Italy 
‘^Amherst Center for Fundamental Interactions and Physics Department, University of Massachusetts, Amherst, MA 01003, USA 
^APC, Univ. Paris Diderot, CNRS/1N2P3, CEA/lrfu, Obs. de Paris, Sorbonne Paris Cite, France 
^Dipartimento di Fisica, Universitd e INFN, Genova 16146, Italy 
^Physics Department, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA 
^INFN Laboratori Nazionali del Gran Sasso, Assergi 67010, Italy 
' Department of Physics, Technische Universita?t Dresden, 01062 Dresden, Germany 
^National Research Nuclear University MEPhl (Moscow Engineering Physics Institute), 115409 Moscow, Russia 
^NRC Kurchatov Institute, Moscow 123182, Russia 
^M. Smoluchowski Institute of Physics, Jagiellonian University, Krakow, 30059, Poland 
Physics and Astronomy Department, University of California Los Angeles (UCLA), Los Angeles, CA 90095, USA 
^Institut fiir Experimentalphysik, Universitdt Hamburg, Germany 
^Physik Department, Technische Universitdt Munchen, Garching 85747, Germany 
^Max-Plank-lnstitut fur Kemphysik, Heidelberg 69029, Germany 
^Dipartimento di Chimica, Universitd e INFN, Perugia 06123, Italy 
^'St. Petersburg Nuclear Physics Institute, Gatchina 188350, Russia 
^Department of Physics, University of Houston, Houston, TX 77204, USA 
^Kiev Institute for Nuclear Research, Kiev 06380, Ukraine 
^Lomonosov Moscow State University Skobeltsyn Institute of Nuclear Physics, Moscow 119234, Russia 
Chemical Engineering Department, Princeton University, Princeton, NJ 08544, USA 
^Gran Sasso Science Institute (INEN), 67100 L’Aquila, Italy 
^Institutfiir Physik, Johannes Gutenberg Universitdt Mainz, 55122 Mainz, Germany 


Abstract 

The Sun is fueled by a series of nuclear reactions that produce the energy that makes it shine. The primary reaction 
is the fusion of two protons into a deuteron, a positron and a neutrino. These neutrinos constitute the vast majority of 
neutrinos reaching Earth, providing us with key information about what goes on at the core of our star. Several ex¬ 
periments have now confirmed the observation of neutrino oscillations by detecting neutrinos from secondary nuclear 
processes in the Sun; this is the first direct spectral measurement of the neutrinos from the keystone proton-proton 
fusion. This observation is a crucial step towards the completion of the spectroscopy of pp-chain neutrinos, as well 
as further validation of the LMA-MSW model of neutrino oscillations. 
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Figure 1: Main nuclear reactions that make up the pp chain (2). The 
neutrinos produced by the chain are highlighted in magenta. Neutri¬ 
nos are named after the parent particles that produce them {pp, pep, 
^Be, ^B). pep, ^Be, and ^B neutrinos have all been previously mea¬ 
sured with Borexino Elllllel, leaving only the pp neutrinos (top left) 
and the extremely faint hep neutrinos (?] (not shown here) undetected 
prior to this work. 


1. Solar neutrinos 

We know that the Sun is fueled by nuclear reac¬ 
tions m. In particular, the “effective reaction” that 
takes place is the conversion of hydrogen (p) into he¬ 
lium (He), 

4p ^'‘He-H 2e+-I-2ve, (1) 

with a net release of kinetic energy and emission of elec¬ 
tron neutrinos Ve- There are two main ways of convert¬ 
ing protons to He nuclei that take place in stars: the pp 
chain and the CNO cycle El. The contribution from 
each of these processes depends on the size, tempera¬ 
ture and age of the star |[T]. The main way of producing 
energy in the Sun and in stars of similar mass and age 
is the pp chain 0. Fig.[T]shows the main reactions that 
are responsible for it. 

The different neutrinos in the chain are often referred 
to by their parent particles. We thus speak of pp neu¬ 
trinos as those produced in the top left of the chain, 
whose parents are simply two protons. Though most of 
the energy and neutrino flux in the Sun comes from the 
pp chain, more massive and hotter stars produce much 
more significant amounts of energy through the CNO 
cycle III. In addition, a small component of the Sun’s 
neutrino flux also comes from the CNO cycle. 


’Now at INFN, 00185 Roma, Italy 

^Now at Georgia Tech, Atlanta, GA 30332, USA 


The challenge in studying neutrinos is intimately con¬ 
nected with the reason why they are of interest in the 
field of solar astrophysics: they can only interact via the 
weak force. Thus, they travel almost undisturbed from 
the core of the Sun to the surface of the Earth, providing 
us with key information about the inner workings of our 
star. 

The energy spectrum of solar neutrinos as predicted 
by the Standard Solar Model 0 (SSM) is shown in 0. 
The end-point energy is E”™ = 420 keV 0. There are 
two currently accepted values for the total pp-neutrino 
flux produced by the Sun, corresponding to the low- 
metalUcity and high-metallicity models 0. The values 
are 5.98 x (1 + 0.006) x 10‘°cm-2 s-’ and 6.03 x (1 + 
0.006) X 10‘°cm“^ s“‘, respectively. 

Various measurements of solar neutrinos 0000 
[la have proven that the Sun is fueled by nuclear reac¬ 
tions and that neutrinos undergo lepton flavor transfor¬ 
mation IT0llII][T2l (also known as oscillations), which 
itself means that at least some of the neutrinos must be 
massive. It also implies that the flux quoted above will 
contain electron-, muon- and tau-neutrinos upon arrival 
at Earth. 

Previous experiments iniiiiiia were able to extract 
the detection rate of pp neutrinos from an integral of the 
solar neutrino spectrum, by subtracting the rates of other 
solar neutrino components cnia. No experiment has 
thus far been able to directly and spectrally detect pp 
neutrinos, that is, to observe interactions of pp neutri¬ 
nos and measure their energies. The main obstacle has 
been the overwhelming background at the low energies 
characteristic of pp neutrinos (see Sec.|^. 

2. Borexino 

Borexino m is a liquid scintillator detector located 
at the Laboratori Nazionali del Gran Sasso, in central 
Italy. It detects neutrinos via their elastic scattering in¬ 
teractions with electrons in the scintillator. Scintillation 
light is collected by photomultipliers (PMTs) surround¬ 
ing the inner detector. Outside the inner detector, a wa¬ 
ter tank as a detector of Cherenkov light, tagging cos- 
mogenic muons with very high efficiency im. 

When combining the pp-neutrino flux provided by 
the Standard Solar Model (Sec. [T]i with the latest values 
of the neutrino oscillation parameters IITSl [191 and with 
information about the Borexino scintillator BOl . we can 
calculate the expected rate of pp-induced electron re¬ 
coils in Borexino to be 131+2 (d-lOOO^' (assuming the 
high-metallicity SSM; the low-metallicity number can 
be obtained by rescaling according to the fluxes). 
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Scintillation light induces signals in ~500 photomul¬ 
tipliers (PMTs) per MeV of deposited energy ED. This 
means that the threshold of the detector can be set as low 
as ~50keV. The end-point of the pp-neutrino energy 
spectrum translates via relativistic kinematics to an end¬ 
point energy of the pp-neutrino-induced electron recoil 
spectrum = 264 keV. Thus, the low energy thresh¬ 
old is especially important for this analysis. 

The pp-neutrino-induced electron-recoil rate is ob¬ 
tained through a ht of the energy distribution of events 
against the known signal and background spectral 
shapes. Our energy estimator is the number of PMTs 
detecting photons within 230 ns after the beginning of 
the physics event. This is converted to energy using 
energy-response and energy-resolution functions tuned 
in previous analyses and re-validated in the current 
study E2l|23]. 

To reduce radioactive background coming from ex¬ 
ternal detector components, we only study events in¬ 
side the centTal fiducial volume (see Sec.|^. The posi¬ 
tions of events are determined by maximizing the like¬ 
lihood of the observed PMT detection time distribu¬ 
tion Eiiiia. The performance of the algorithm has 
been tested with calibration sources EH; further de¬ 
tails can be found in Il22l . At the energies relevant for 
the pp analysis, the position reconstruction code has an 
uncertainty of 20 cm. This was estimated by compar¬ 
ing the reconstructed and nominal positions of radioac¬ 
tive sources placed inside Borexino during a calibration 
campaign ll^l^. 

The main backgrounds in this analysis are in¬ 
trinsic to the organic liquid scintillator, and its pile- 
up. Other backgrounds include ^'‘’Po, *^Kr and 

^'"^Pb. The rate of /3-decay in Borexino is ~ 40 (s ■ 
loot) * II 20 I . four orders of magnitude larger than the 
expected pp-induced recoil rate (above). The strategy 
for assessing the real background rates is discussed in 
Sec. [31 

3. pp analysis 

Data for the present analysis was collected between 
January 2012 and June 2013, containing the hrst 408 
live days of data of Borexino Phase-II. These data are 
characterized by reduced levels of some of the most 
prominent radiogenic backgrounds, most notably ^^Kr, 
^'°Bi and ^'°Po, thanks to a series of successful purihca- 
tion campaigns Il22l . Basic selection criteria have been 
described previously in EH, in the context of Borexino 
Phase-I. In particular, the criteria in the current analysis 
are: 

(1) No coincidence with muon events; 



Figure 2: Spectra of first and second events in trigger windows ac¬ 
quired during the live time considered in the pp analysis, with the 
cuts outlined in the text. Second clusters are not subject to the trigger 
threshold characteristic of Borexino Ga, and are thus useful for the 
independent estimate of the rate. 


(2) Event triggered by a scintillation event in the inner 
detector (Sec.j^; 

(3) No electronic noise from PMTs or electronics racks; 

(4) No coincidence with isotopic decays in the radon 
branch of the uranium series; 

(5) Position reconstruction within the hducial volume; 
for this analysis, the hducial volume is a truncated 
sphere with R - 3.021 m, |z| = 1.67 m; the hducial vol¬ 
ume is 86 m^, or 75.51. 

As explained in Sec.[^ most of the remaining back¬ 
ground is ;0-decays. To estimate its rate indepen¬ 
dently from the main analysis, we looked at a sample of 
data in which the event causing the trigger is followed 
by a second event within the same trigger window (16- 
ps long, while physical events last ~ 100-500 ns). These 
second events are not subject to the trigger threshold, so 
they register down to much lower energies. The spectra 
of hrst and second events are shown in Fig.|^ By htting 
this spectrum against the known spectral shape Il27l . 
we obtain a rate in Borexino of 40+1 Bq per 1001. 
The error accounts for statistical and systematical un¬ 
certainties, most notably due to uncertainty in the 
spectral shape Il271l28lI29ll and ht conditions. 

Due to its high rate, naturally generates a con¬ 
siderable amount of pile-up: two physical events occur 
so close to each other in time that the detector fails to 
resolve them. With the above measurement of the 
rate, the expected rate of pile-up is ~100 per 

day per 1001. This is comparable to the expected pp 
neutrino interaction rate (Sec. j^. In addition, the end¬ 
point of the pile-up spectrum is 312keV (twice 

the spectrum end-point EOj), also close to the end- 
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point of the pp-induced recoil spectrum, This un¬ 
derlines the importance of a careful assessment of pile- 
up events, both and other pile-up. 

The pile-up component can be determined using an 
independent, data-driven method called synthetic pile- 
up. Real triggered events with no cuts are artificially 
overlapped with random data obtained from the ends of 
real trigger windows, uncorrelated with the triggering 
event. By overlapping four random data samples with 
each physical trigger, we obtain a sample of synthetic 
events that has four times the exposure of the real data 
set used. The synthetic events are then reconstructed 
using the same software used for real events, and se¬ 
lected with the same criteria. Using this method, we ob¬ 
tained the true rate and spectral shape of pile-up in our 
detector. This method naturally accounts for the effect 
of fiducialization on pile-up events. 

With the and pile-up rates fixed, we obtained the 
pp rate by fitting the data spectrum in a window of 165- 
590 keV against the known spectral shapes. The fit was 
performed with a tool named spectral-fitter, previously 
used in ^Be jS) and pep @ analyses. The tool was im¬ 
proved to work on any Borexino analysis. In particu¬ 
lar, we included a treatment of the scintillator energy 
response to mono-energetic electrons at high statistics; 
a modified treatment of the energy resolution at low en¬ 
ergy; and the introduction of synthetic pile-up. 

In addition to pp, and pile-up, other species in¬ 
cluded in the fit are other solar neutrinos (^Be, CNO and 
pep), *^Kr, ^'°Bi, ^'°Po and The ^Be rate was 

fixed to the rate previously obtained by Borexino 0; 
CNO and pep rates were fixed to the values predicted 
by the high-metallicity Standard Solar Model; remain¬ 
ing background rates are left free in the fit, except for 
^'"^Pb, which is fixed to the value measured by looking 
at ^''^Bi-^^'^Po coincidences Il22l . The scintillator light 
yield and two energy resolution parameters are free to 
vary in the fit. 

The energy spectrum after cuts, together with the best 
fit, is shown on Fig. The uncertainties shown in the 
figure are statistical only. 

We considered an alternative method to account for 
pile-up, known as convolution method. Random 16-jUS- 
long samples of data are regularly collected by Borex¬ 
ino ED. These data are sliced into time windows of 
the same duration as the window used for our energy 
estimator (Sec.|^. By computing the energy estimator 
inside each of these small windows, we generate the dis¬ 
tribution shown in Fig. This distribution is the spec¬ 
trum collected by Borexino without a trigger threshold, 
and it includes electronic noise, as well as all physi¬ 
cal events generated by signal and background species. 



Energy (keV) 


Figure 3: a. Fit of the data in our energy region of interest to extract 
the value of the pp neutrino interaction rate. The ^Be rate displayed 
corresponds only to the higher-energy branch O. Rates are displayed 
in (d-100t)“^ except for that of which is in (s-100t)“^ Uncer¬ 
tainties are statistical only. Species labelled “constrained” are free to 
vary, but a penalty is applied to the likelihood if the value chosen by 
the fitter differs from the expected value (23l. b. Residuals of the fit. 
Systematic effects are discussed in the text. 



0 50 100 150 200 250 


Energy estimator: number of hit PMTs 

Figure 4: Spectrum obtained by slicing random triggers into smaller 
time windows with the same duration as the windows used in regu¬ 
lar data for the computation of our energy estimator, as explained in 
Sec.[^ This is the spectrum obtained by Borexino without a trigger 
threshold. It is used in the current analysis for an alternative estimate 
of the pile-up rate. 
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Figure 5: pp neutrino interaction rate as obtained from many fits 
where fit conditions are varied to account for the systematic effects 
outlined in the text. The mean value and root mean square of this 
distribution are our central value and systematic uncertainty, respec¬ 
tively. 


Pile-up can be generated by convolving each spectral 
component (neutrinos and backgrounds) with this spec¬ 
trum. 

Fit conditions were varied to estimate the systematic 
uncertainty. Effects considered include variations in; 

(1) The pile-up evaluation method (synthetic and convo¬ 
lution)', 

(2) The length of the time window in which we count 
PMTs hit (Sec.|5; 

(3) The definition of the fiducial volume; 

(4) The expected value of the CNO neutrino detection 
rate; 

(5) The fit start and end points. 

Fits were performed using many different combinations 
of all these effects. The histogram of the rates obtained 
in all these conditions is shown in Fig.|^ The root mean 
square of the distribution is 9.8 (d-1001) '. In addition, 
the resolution of the position reconstruction algorithm 
(20 cm; Sec. ^ induces a 2% uncertainty in the deter¬ 
mination of the fiducial volume. These uncertainties are 
added in quadrature to obtain the total systematic un¬ 
certainty of our measurement. The final pp neutrino 
interaction rate is m 

Rpp = 144+ 13(stat) + 10.(syst)(d- 100 0^' (2) 

We also studied possible variations in the result caused 
by; 

(1) Uncertainties in the spectral shapes of 127] [28] 
EH and^iOBi ||22l[32l[33l ; 

(2) Variations of the expected pep rate within the mea¬ 
sured and theoretical uncertainties; 


(3) Leaving the rate free in the fit; 

(4) Varying the veto time after muons cross the detector. 
No significant variations were observed. 

We also studied as a potential source of background 
the j8-emitter *^Rb, with a Q-value of 283.3 keV Il30l . 
very similar to the end-point of the pp-induced elec¬ 
tron recoil spectrum (Sec. |^. Using the known 
crustal abundances 1341 . isotopic abundances and half- 
lives l30ll of *^Rb and we can calculate the ratio 
between their natural activies to be 1 to 4, respectively. 
Under the assumption that the relative abundances of 
^^Rb and in Borexino are equal to their relative 
crustal abundances, and knowing that the rate of in 
Borexino is < 0.4(d-100t) * fU, we can constrain the 
^^Rbrateto <0.1 (d-1001)-'. 

4. Interpretation 

The electron neutrino survival probability. Pee, is the 
probability for a solar neutrino produced as electron- 
type to arrive at Earth as an electron-neutrino. It can 
be computed using the relations in l35l [3^ . the most 
up-to-date oscillation parameters and interaction cross- 
sections from lT8l[T9ll37ll and the electron number den¬ 
sity in Borexino, (3.307±0.003)xl0^' (lOOt)-' fS). Us¬ 
ing that number and our measured pp neutrino inter¬ 
action rate, we can calculate the pp neutrino flux to be 
(6.6 ±0.7)x 10'° cm“^ s“', consistent with the SSM pre¬ 
diction (see Sec.[T]l. The uncertainty of the present mea¬ 
surement does not allow us to distinguish between the 
high- and low-metallicity solar models. 

Conversely, if we use the pp neutrino flux predicted 
by the SSM as an input, we can calculate the survival 
probability Pee - 0.64 + 0.12, consistent with expec¬ 
tation from the LMA-MSW solution ll38]l . as shown in 

FigJU 

Finally, the fact that the present measurement of the 
pp neutrino flux is consistent with expectation justifies 
the assumption made by the SSM that the photon and 
neutrino luminosities from the Sun can be correlated; 
i.e., that the Sun has been stable for the past 10^ years. 
This invalidates a previous hypothesis that changes in 
the solar core over these time scales could be responsi¬ 
ble for terrestrial global epochs Il39l . 
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